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Abstract. In single-field infiationary models with a low sound speed, the orthogonal 
shape of the primordial bispectrum arises due to partial cancellations between equilateral- 
type shapes. This fact allows for a speed of sound Cg as low as about 0.01, which is actually 
weakly preferred by WMAP data. For such values, the trispectrum, scaling like 1/cf, is 
of order 10^ and is therefore comparable to, and greater than, the la observational bound 
~ (—3.11 ± 7.5) X 10^. Hence, the trispectrum is already constraining infiationary 
mechanisms candidates for generating an orthogonal bispectrum at the level hinted in 
WMAP data. If this signal persists in imminent Planck data, most of the parameter space 
of the simplest effective field theory of infiation will be under observational pressure, while 
a dedicated analysis will be needed for the substantial fraction of parameter space where 
we show that a qualitatively new, orthogonal, trispectrum naturally arises. 
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1 Introduction 

The deviation from perfect Gaussian statistics of the primordial curvature perturbation 
( enables one to discriminate amongst the candidate physical mechanisms that produced 
the seed primordial fluctuations (see for instance [1-5] for recent reviews). In this respect, 
it is fair to say that, despite significant efforts, the trispectrum has received considerably 
less attention than the bispectrum, in particular concerning data analysis. To our knowl- 
edge, three types of constraints on well motivated primordial trispectra are now available^: 
constraints on Tjy^ and g^i [6, 7] - which set the amplitude of the two different trispectra 
generated classically on super-Hubble scales in early-universe models with multiple de- 
grees of freedom - and constraints on [6], setting the amplitude of a representative 
'equilateral-type' trispectrum, generated by quantum interactions around Hubble crossing. 
The latter parameter is the equivalent for the trispectrum of f^j^ for the bispectrum, while 
the former are the counterparts of Another important type of primordial bispectrum 
constrained by data, orthogonal non-Gaussianities [8], has up to now no counterpart at 
the level of the trispectrum. 

The purpose of the present paper is two-fold. Our first aim is to show that an 
'orthogonal-type' trispectrum naturally arises in a significant fraction of parameter space 
in the simplest theoretical context (to be more accurate, this is a one-parameter family of 

^Reference [6] constrains as well the non-primordial trispectrum signal induced by cosmic strings, as 
well as the primordial constant model, which provides a useful benchmark, but which has no physical 
motivation yet. 
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trispectra as we will see). Our second, related, aim, is to point out the use of the trispec- 
trum as a useful diagnostic for the appearance of a large bispectrum of the orthogonal 
type. 

In particular, the final WMAP data [9] indicated a 2.45cr hint of orthogonal non- 
Gaussianities: —445 < f^j^ < —45 (95% CL) (while showing no evidence of equilateral 
non-Gaussianities: —221 < < 323 (95% CL)). When interpreted in terms of the 
effective field theory of inflation [8, 10], the orthogonal shape arises from partial cancella- 
tions between otherwise equilateral shapes^ (see [12, 13] for the first concrete realization 
of this mechanism), leading to a smaller amplitude than the general estimate Jni ~ l/c^, 
namely f^f^ = C(^). This fact allows for a speed of sound Cg as low as about 0.01, 
which is actually (weakly) preferred by current data [9] (see Fig. 1). For such values of c^, 
and unless inflation occurred in the region of parameter space where similar partial can- 
cellations that leads to the orthogonal bispectrum arises at the level of the trispectrum, 
the amplitude of the trispectrum, scaling as 1/c^, is of order 10® and hence is already 
comparable to (and greater than) current constraints t']^j^ = (—3.11 ± 7.5) x 10^ obtained 
with WMAP data [6] . We will make this more quantitative in the body of this paper, but 
our message is simple: the trispectrum is already constraining inflationary mechanisms 
candidates for generating an orthogonal bispectrum at the level hinted in WMAP data. If 
this signal persists in imminent Planck data, most of the parameter space of the simplest 
effective field theory of inflation will be under strain, while a dedicated analysis of our 
orthogonal trispectrum signal will be needed for the remaining one. 

The plan of our paper is as follows. In section 2, we use the effective field theory 
of inflation at the single-derivative level to parametrize the cubic and quartic action for 
fluctuations in low sound speed models. We then give the expression of the trispectrum 
generated in these models, using the results of Chen et al [14] in the setup of k-inflation 
[15, 16], which is computationally equivalent. The overall amplitude of this trispectrum 
is fixed by the speed of sound c^, while its shape depends on two-parameters: A, which 
determines the shape of the bispectrum, and B, which is unrelated in general. Section 3 is 
then dedicated to the study of the resulting 2-parameter family of shapes of trispectra in 
the region of parameter space 3.1 < A < 4.2 where the orthogonal bispectrum is generated 
and where values of the speed of sound as low as (9(0.01) are allowed. We define the region 
of parameter space where the trispectrum can be well represented by the 'equilateral' one 
and is thus already constrained by present data, and represent the new shape that arises 
in the complementary region. We conclude in section 4. 

■^Recently, a bispectrum with a significant overlap with the orthogonal shape was shown to arise in a 
different context [11]. 
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2 The trispectrum from low sound speed models 



In this section, we give the expression of the leading-order trispectrum generated in the 
simplest set-up of the effective field theory of inflation, which is computationally equivalent 
to k-inflation. Readers who are familiar with both the effective field theory of inflation and 
the trispectrum generated in low sound speed models can skip this section and proceed 
directly to section 3. 

2.1 The effective field theory of inflation up to quartic order 

In this subsection, we briefly review the effective fleld theory of inflation developed in 
[10], or to be more accurate the effective fleld theory of fluctuations generated in single- 
clock inflation. In such models with only one non-gravitational degree of freedom, it is 
always possible to choose a slicing such that surfaces of constant t coincide with surfaces 
where the 'clock' is unperturbed, i.e. such that 6(j){t,x) = in the case of a scalar 
fleld clock. No explicit scalar fluctuations appear in this unitary gauge in which time 
diffeomorphisms have been flxed. The most generic effective action in this gauge can thus 
be built by allowing only metric operators invariant under the unbroken time- dependent 
spatial reparametrizations. It can then be shown that considering fluctuations around a 
spatially flat FLRW background amounts to studying the following action [10] 



S = d Xa/— ^ 



where H is the Hubble parameter, 5(7°° = (?°° + 1, SK^^, (respectively SR^yp^) is the 
fluctuation of the extrinsic curvature of constant time surfaces (respectively of the 4- 
dimensional Riemann tensor) and where F starts quadratic in its arguments 5g^^, ^K^^, 
and 5R^ypa. The simplest effective fleld theory at lowest order in derivatives corresponds 
then to allowing operators involving powers of 5g^^ only, namely, up to quartic order in 
fluctuations, 

F = ^M,itn5g'y + ^M,{tn5g'y + ^M.mSg'y . (2.1) 

The true effectiveness of this approach relies on the gravitational analogue of the equiva- 
lence theorem for the longitudinal components of a massive gauge boson [17]: the scalar 
degree of freedom can be explicitly reintroduced in the Stiickelberg trick, thus restoring 
full time-diffeomorphism invariance, and it decouples from the gravitational sector at high 
enough energies, allowing to neglect the complications of the mixing with gravity. In this 
decoupling regime, the effect of the Stiickelberg time diffeomorphism t — t- t + 7[{x) on Sg'^^ 
is simply 

5^00 _ ^2 ^ , (2.2) 
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while one can neglect the terms introduced by the time dependence of the M„(t) at leading 
order in a slow-varying approximation, or equivalently by assuming that vr enjoys an 
approximate shift-symmetry. The resulting Lagrangian reads, up to quartic order in vr: 



^1 



DL 



d Xy 



-9 



+ 



4 



i2\2 



/ ..3 , O.2.;. ^2^ , 2M4^.4 



(2.3) 



where (^^vr)^ = — tt^ + (5j7r)^/a^ is evaluated on the background metric and vr is related to 
the curvature perturbation by the simple relation ( = —Htt at linear order and at leading 
order in a slow-varying approximation. A non-zero M2 introduces both a reduced sound 
speed Cs, such that 

1 



2M| 



(2.4) 



and large cubic and quartic interactions. By introducing A/cj. = —1 
the WMAP notation in [9] , and 
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the decoupling Lagrangian can be cast in the form 

,(9.vr)2 
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(2.6) 



where we have kept leading-order terms in the interesting limit <^ 1 of a large bispec- 
trum. As explained in [8], A of order one is technically natural from the effective field 
point of view as the operators in tt^ and nldiir)"^ then introduce the same strong coupling 
scale. The same reasoning shows that B of order one is technically natural as well. Note 
also that, as the effective field theory of inflation with operators involving only powers 
of Sg^^ is computationally equivalent to k-inflation, it should not come as a surprise that 
the Lagrangian (2.6) can be identified with the one in [14, 18] (see also [19]), with the 
correspondence, at leading order in 1/c^: 



A A 

^ ^ -2- 



B a 



9: 



(2.7) 



Note eventually that DBI inflation [20, 21] simply corresponds to A 
parametrization. 



-B 
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2.2 The trispectrum 

Using the correspondence (2.7), the primordial trispectrum generated from the inflationary 
fluctuations Lagrangian (2.6) can be simply read off from the equivalent result Eq. (3.32) 
in [14], namely 



4 1 

(C(ki)C(k2)C(k3)C(k4))c = {2Txf Vl5{Y, ki) n h, h2, h,) (2.8) 

i 

where kij = \ki + kj\, 



1.3 

i=l 



r(A, 5) = i (^T,i - ^T,2 + r,3 - BT,, ) , (2.9) 



{ki + k2 + h + hy' ^ ^ 

and where the explicit (lengthy) expressions of Tsi,Ts2,Ts2, can be found in [14]. Note 
that although the operators in ((SjTr)^)^ and Tr^{di7iY in Eq. (2.6) are of the same order as 
the one in tt^, they cancel in the quartic Hamiltonian at leading order in the small sound 
speed limit, leaving only the scalar exchange contributions T^j and the contact-interaction 
trispectrum generated from the operator in tt^. 

A detailed analysis of the shapes of the four constituent trispectra in (2.9) was per- 
formed in [14], to which we refer the reader for more details (see also [19]). Overall, 
they reached the conclusion that they all share very similar properties. In particular, as 
expected from trispectra of quantum origin generated around the time of Hubble cross- 
ing, they are the largest for the conflgurations where both the external and internal 
momenta are of similar magnitude, i.e. near the regular tetrahedron (RT) limit where 
ki = k2 = ks = = ki2 = ki4. This similarity was used in [6] in which observational 
constraints on the simple and representative 'equilateral-type' trispectrum Td were de- 
rived (the only observational constraint on a trispectrum from quantum origin to date). 
However, we are interested here in the possible cancellations between these overall simi- 
lar shapes, as we vary the parameters {A,B), which could result in a trispectrum poorly 
correlated with Td, and to which current constraints would hence be blind. As we have ex- 
plained in the introduction, of particular interest is the region (weakly) favored by WMAP 
nine-year data 3.1 < A < 4.2 where partial cancellations between the operators 7r{di7r)'^ 
and TT^ in (2.6) leads to a primordial bispectrum correlated with the orthogonal template 
at more than 80 % [8], and in which a low sound speed of order 0.01 is allowed [9] (see 
Fig. 1). We present the results of such a study in the following section. 
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Figure 1. WMAP nine-year constraints: la, 2a and So" confidence regions derived from the 
bispectrum on the sound speed Cg and interaction coefficient A in Eq. (2.6). Figure taken from 



3 'Equilateral' and orthogonal trispectra. 
3.1 Results 

As we have just explained in the preceding section, we would like to assess where in 
parameter space the overall trispectrum (2.9) can/cannot be well represented by the simple 
'equilateral-type' trispectrum Td constrained by data. Given the expression (2.9), it is 
clear that, for any A, the trispectrum interpolates between highly correlated and highly 
anti-correlated with as we vary B from large negative to large positive values. A 
qualitatively new shape should hence arise in the neighborhood of a particular value of B. 
The only questions are then: around which value? How narrow is this region? What the 
new shape looks like? And how much does the latter vary with AP. in particular in the 
region 3.1 < A < 4.2. 

To answer these questions, one should ultimately resort to the scalar product defined 
by the estimator used to constrain the primordial trispectra, for instance the correlator 



[9]. 



[22, 23] 




(3.1) 
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between angular trispectra when using the CMB as the observational probe. Implement- 
ing such a correlator is however beyond the scope of this paper. We rather sticked to 
representing 2-dimensional slices of our trispectra in different representative limits, as was 
originally done in [14]. Note that we actually did make use of quantitative correlators, in 
Fourier-space, using reduced trispectra [23, 24] or full trispectra [25]. However, we found 
these correlators to be somewhat misleading for our purpose, as they indicated widely 
different regions of parameter space as giving rise to an 'orthogonal' trispectrum, in the 
sense of a trispectrum with a small correlation with Td. Moreover, when visually rep- 
resenting these candidate orthogonal trispectra, they appeared almost indistinguishable 
from Tel. On the contrary, we believe that our procedure is trustworthy and sufficient to 
define where represents well or not the overal trispectrum (2.9), and, as we will explain 
below, we have quantitative arguments that support our findings. 

To summarize the latter, we find that, for the three representative values A = 
{3.2,3.6,4} in the interesting range 3.1 < A < 4.2, the total trispectrum (2.9) can be 
well represented by T^i only for 

• < 5 and > 11 (A = 3.2) (3.2) 

• B<5 and B>U{A = 3.6) (3.3) 

• 5 < 10 and > 19 (A = 4) (3.4) 

while a qualitatively different shape arises in the complementary region, centered around 
the values 

• ~ 8.5 (A = 3.2) (3.5) 

• ~ 11.5 (A = 3.6) (3.6) 

• ~ 14.5 {A = 4) (3.7) 

Moreover, we find that the shapes of these various orthogonal trispectra depend very 
weakly on A and can thus be well represented by (up to an overall multiplicative factor) 

Torth = 3.2 T,i - 1.8T,2 + T,3 - 11.5 T,i , (3.8) 

corresponding to the values A = 3.6, B = 11.5. 



Quantitatively, a simple measure of the amplitude of the trispectrum is given by the 
parameter t^L defined such that 

r^T{ki, k2, h, k4, ku, ku) Inl , (3.9) 

K'^ limit 
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where RT stands for the regular tetrahedron hmit ^1 = ^2 = ^3 = ^4 = ki2 = ki^ = k. 
Apphed to (2.9), this gives 

t^L = \ (0.062A^ - 0.210A + 0.305 - 0.0355) . (3.10) 

The observational constraint relevant for our purpose in [6] is derived by assuming that T 
in Eq. (2.8) can be well approximated by in (2.10), i. e. it puts a bound on the parameter 
defined such that T = f^i/ttqiiTd) x T^. It is therefore applicable only in the ranges 
defined in Eqs. (3.2)-(3.4) where the trispectrum becomes essentially proportional to Td. 
To make the link with the observational constraints on the bispectrum, let us apply this 
to the value of the speed of sound = 0.013^ such that f°^f = -245 (the central WMAP 
estimate) for the representative value A = 3.6. One then finds 

tjvL = 1.17 X 10^ X (1 - 0.0985) (3.11) 

which is comparable to (and actually greater than) the la constraint [6] 

= (-3.11 ±7.5) X 10^ (3.12) 

Similar numbers are of course found in the range 3.1 ^ A < 4.2: Inl = 1.62 x 10^ x (1 — 
0.0985) for A = 3.2 and tjvL = 8.80 x 10^ x (1 - 0.0985) for A = 4, under the same 
hypotheses. There is currently no point in performing a more detailed statistical analysis 
given the weak 2.4:5a hint of an orthogonal bispectrum, but we believe our message is 
clear: the trispectrum is already constraining candidate low sound speed infiationary 
models generating an orthogonal bispectrum at the level suggested in WMAP data. For 
instance, the result Eq. (3.11) is 2.13a (respectively 3.50ct) away from the central value 
Eq. (3.12) for 5 = —1 (respectively 5 = —10). Even a decrease in the error bar by a 
factor of a few, as expected from Planck [26], could hence constrain these scenarios in a 
statistically significant way, provided of course that the orthogonal bispectrum signal is 
confirmed. 

Note once again that Eq. (3.11) can be meaningfully compared to the observational 
constraint (3.12) only for 5 < 5 and 5 > 14. In the intermediate range, a shape quali- 
tatively different from Td arises to which current constraints are mostly blind. Relatedly, 
the fact that t^i in Eq. (3.11) vanishes for 5 = 10.20, so in the middle of this intermediate 
range, and close to the value 5 = 11.5 at which we defined our representative orthogonal 
trispectrum, justifies a posteriori our procedure and our findings. Indeed, t^i measures 
the amplitude of equilateral- type shapes, which peak around the regular tetrahedron limit. 
The fact that it vanishes does not indicate that no significant non-Gaussianities are gen- 
erated, but rather that they can not be faithfully represented by the 'equilateral' ansatz 

''This value of Cg is found using Eq. 57 in Rcf. [9]. 
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Tel in Eq. (2.10) The estimator t^vL in Eq. (3.10) vanishes as well at -B = 7.76 for A = 3.2 
and at -B = 13.21 for A = 4, so again in the intermediate ranges that we defined and close 
to the values B ~ 8.5 and B ~ 14.5 at which defined the appearance of the orthogonal 
trispectrum. 

Eventually, let us stress that the origin of our representative orthogonal trispectrum 
Eq. (3.8) is not an ad-hoc orthogonalization procedure: by subtracting out the similari- 
ties between Td, T^i, Ts2 and T^s, one can indeed imagine the construction of a basis of 
the vector space spanned by these four trispectra constituted by Td and three qualita- 
tively different and mutually orthogonal trispectra. While mathematically correct, this 
procedure would be somewhat artificial. On the contrary, we have shown the natural ap- 
pearance of trispectra qualitatively different from the equilateral one Td, and which can 
be represented by the template Eq. (3.8), in a substantial fraction of parameter space in 
the simplest theoretical context. 

3.2 Shapes of trispectra 

In this subsection, we show plots of Td and of T(3.6, B) supporting the results stated 
above, namely that the latter can be represented to a good approximation by Td for 
B < 5 and B > 14, and that a trispectrum qualitatively different from Td arises near 
B = 11.5. We leave other plots, in particular for other values of A, to the appendix A. 

Our scale invariant trispectra are in general functions of 5 variables, and in what 
follows, we follow Ref. [14] and represent 2-dimensional slices of them in different repre- 
sentative limits. A slight difference is that we chose to plot the scale-independent quantities 
T = T/ {kik2kskiY^'^ rather than T itself (we have checked that similar conclusions are 
reached by using the two sets). The four limits we consider are: 

• The specialized planar limit, in which ki = ks = ku, and the tetrahedron reduces 
to a planar quadrangle with [14] 



Shapes are then represented as functions of /c2/^i and k^/ki. 

• Near the double-squeezed limit: k^ = k^ = ki2 and the tetrahedron is a planar 
quadrangle with [14] 




(3.13) 



, (_1.2 _i_ 1.2 I t.2\ _ i,2 1,2 I 1.2 1.2 i U2 t.2 i t.2 t.2 _ 1.2 l2 _ iL4 , l2 l2 

V2k4 



(3.14) 
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where kgi and are defined as 



fc^i = 2a/ {kik4^ + ki ■ k4)(/ci/c4 - ki ■ , 



ks2 = 2a/ (/i;3A;4 + ks ■ k4)(A;3A;4 - ks ■ k4) . 



(3.15) 



Shapes are then represented as functions of k^/ki and ki^/ki. 

• The folded limit: ki2 = 0, hence ki = k2 and ^3 = /C4. Shapes are then represented 
as functions of k^/ki and ki^/ki, and we assumed A;4 < ki without loss of generality. 

• The equilateral limit: ki = k2 = k^ = k^. Shapes are then represented as functions 
of ki2/ki and ku/ki. 

In Figs. 2, 3 and 4, we represent the shape functions of (top left), T(3.6,11.5) 
(top right), T(3.6,5) (bottom left)) and — T(3.6, 14) (bottom right), in the specialized 
planar limit, near the double-squeezed limit and in the folded limit respectively. The 
shape functions are left white when the momenta do not form a tetrahedron. From these 
plots, it is evident that T(3.6, 5) (respectively T(3.6, 14)) is well correlated (respectively 
anti-correlated) with T^, while T(3.6, 11.5) is qualitatively different from it. Note that 
the scales are different in each plot, and that T(3.6, 11.5) has both positive and negative 
values. In these three limits, the tetrahedron reduces to a planar quadrangle, which is the 
configuration probed by the CMB (see [23]). In Fig. 5, the non-planar, equilateral, limit, 
is also displayed. T^i is constant in this limit, so that its effect in Eq. (2.9) is simply an 
overall shift in amplitude. For this reason, we represent only T^i (left) and T(3.6,11.5) 



Eventually, note that for A < 3.1 and A > 4.2, the bispectrum is of equilateral type 
and observationally allowed values of the speed of sound Cg are greater than in the range 
3-1 ^ A < 4.2, of order > 0.1 [9] (see Fig. 1). The trispectrum (2.9) is then of maximum 
amplitude C(10'*) and, for such values, is not expected to be efficiently constrained by CMB 
data [27]. It is nonetheless interesting to investigate the appearance of orthogonal- type 
trispectra in this region of parameter space. We have not performed an exhaustive study 
of this, but the reasonings made in subsection 3.1 show that a trispectrum qualitatively 
different from T^i arises around a certain value of B for any A, and that the cancellation 
of the estimator ^atl in Eq. (3.10) is a good estimate of where this arises, leading to the 
one-parameter family of approximate orthogonal trispectra 



We have visually checked that this indeed the case, and moreover that outside the range 
3-1 ^ ^ ^ 4.2, these trispectra can differ substantially from our representative orthogonal 
trispectrum Eq. (3.8). 



(right). 




(3.16) 
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Figure 2. In this group of figures, we consider the speciahzed planar hmit with ki = = /ci4, 
and plot fci (top left), f'(3.6, 11.5) (top right), r(3.6, 5) (bottom left) and - f{?,.Q, 14) (bottom 
right) as functions of ^2/^1 and k^/ki. 

4 Conclusions 

Amongst the three template bispectra constrained by the WMAP team, the orthogonal 
shape has been much less studied than local or equilateral non-Gaussianities. Although 
not statistically significant, the 2.4:5a hint of orthogonal bispectrum indicated by the 
final data [9] hence deserves close attention. Observational consistency checks have been 
performed in [9], showing no obvious source of systematic contamination. The purpose 
of this paper was to show that the trispectrum is already providing another important 
consistency check, and that it actually puts under strain some of the simplest inflationary 
mechanisms generating such a signal. The reason is simple: when interpreted in the 
theoretical framework of the effective field theory of single-clock inflation, the WMAP 
signal favors a value of the sound speed of inflaton fluctuations of order (9(0.01). For 
such low values, the trispectrum, scaling like 1/c^, is generically of order 10^, and hence 
is comparable to, and greater than, the current la observational bound t^^j^ = (—3.11 ± 
7.5) X 10^ [6]. 

This limit is derived under the assumption that the trispectrum generated in low 
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in 



Figure 3. In this group of figures, we look at the shapes near the double squeezed limit: we 
consider the case where = = ki2 and the tetrahedron is a planar quadrangle. We plot 
Tci (top left), f (3.6, 11.5) (top right), f (3.6, 5) (bottom left) and -f'(3.6, 14) (bottom right) as 
functions of k/^/ki and ku/ki. 

sound speed models can be well approximated by the simple 'equilateral-type' shape gen- 
erated by the operator tt^. While this is true in a large fraction of parameter space, we have 
shown that there exists an non-negligible fraction of it in which qualitatively new shapes 
naturally arise, to which current constraints are almost blind. In the interesting region in 
which an orthogonal bispectrum is generated, an essentially unique orthogonal trispectrum 
Eq. (3.8) can be generated, and it would be interesting to undertake an analysis of CMB 
data in search for such a signal. 
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Figure 4. In this group of figures, we consider the folded limit ki2 = 0, and plot Td (top left), 
^(S. 6, 11.5) (top right), T(3.6, 5) (bottom left) and — T(3.6, 14) (bottom right) as functions of 
/C4//C1 and ku/ki. 




Figure 5. In this group of figures, we consider the equilateral limit ki = k2 = = k^, and 
plot Tci (left) and T"(3.6, 11.5) (right) as functions of ^12/^1 and ku/ki. 
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A Comparison between orthogonal trispectra 

We have explained in subsection 3.1 that the trispectrum Eq. (3.8) represents well the 
orthogonal trispectra Eqs. (3.5)-(3.7) arising in the range S.l < A < 4.2. In this appendix, 
we demonstrate this by showing, in Figs. 6, 7, 8 and 9, plots of T(3.2,8.5), T(3.6,11.5) 
and T(4, 14.5) in the same limits as for Figs. 2 to 5. For comparison, we also represent 
Tgi, an 'equilateral- type' shape different from T^, which we have represented in the other 
figures. The similarity between 7^(3.2,8.5), T(3.6, 11.5) and 7^(4, 14.5) is striking. 




Figure 6. In this group of figures, we consider the specialized planar limit with ki = = A;i4, 
and plot fsi (top left), r(3.6, 11.5) (top right), 'f(3.2,8.5) (bottom left) and r(4, 14.5) (bottom 
right) as functions of k2/ki and k^/ki. 
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